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Nuclear masses are calculated using the modified Bethe-Weizsacker mass formula in which the iso- 
tonic shifts have been incorporated. The results are compared with the improved liquid drop model 
with isotonic shift. Mass excesses predicted by this method compares well with the microscopic- 
macroscopic model while being much more simple. The neutron and proton drip lines have been 
predicted using this modified Bethe-Weizsacker mass formula with isotonic shifts. 
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Nuclear reactions using radioactive ion beams have been studied extensively to allow measurements of the masses, 
half-lives, radii and other properties of unstable nuclei. These results reveal some interesting characteristics of unstable 
nuclei, such as a halo and a skin [1]. The disappearance of some traditional magic numbers and extra stability for 
some neutron numbers [2] have also been observed. The Bethe-Weizsacker (BW) mass formula designed to reproduce 
the gross features of nuclear binding energies for medium and heavy nuclei fails for light nuclei, specially away from 
the line of stability [3]. The unusual stability of nuclei with preferred nucleon numbers, commonly referred to as 
magic numbers, can be clearly delineated by comparing the experimental binding energies with predictions of the 
Bethe-Weizsacker mass formula. However, this comparison does not indicate the recently observed features like the 
disappearance of some traditional magic numbers and extra stability for some newly observed nuclei. To alleviate these 
problems, a generalized Bethe-Weizsacker mass formula was predicted [4] which explained the gross features of the 
shapes of the binding energy versus neutron number curves of all the elements from Li to Bi. Both Bethe-Weizsacker 
and the modified Bethe-Weizsacker (BWM) mass formulae have no shell effect incorporated and thus when the shell 
effect in a nucleus quenches, the theoretical mass comes closer to the experimental ones. A comparison with the 
experimental data and the BWM delineates the positions of all old and the new magic numbers. For heavier nuclei 
this BWM mass formula approaches the old BW mass formula. 

In this work we improve the mass predictions of BWM by introducing the isotonic mass shifts [5]. We find that 
it greatly improves the agreement between the experimental masses and the calculated ones and thus predicts the 
driplines more accurately than calculated earlier [6] with BWM alone. We adopt a procedure in which we compute 
the average shift of the BWM formula from the experimental values. The shift, called the isotonic shift A„, is then 
subtracted from the prediction of the BWM formula [5], as explained later. 

In the liquid drop description of a nucleus of mass number A and atomic number Z, the binding energy is given by, 

B(A, Z) = a v A - a s A 2 / 3 - a c Z(Z - l)/^ 1 / 3 - a sym (A - 2Z) 2 /A + 5, (1) 

where 

a v = 15.78 MeV, a s = 18.34 MeV, a c = 0.71 MeV, a sym = 23.21 MeV and a p = 12 MeV, (2) 
and the pairing term, 

S = a v A~ x l 2 for even N — even Z, 
= -a p A- 1/2 for odd N - odd Z, 
= for odd A, 

(3) 
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The above formula was prescribed by Bethe-Weizsacker (BW). With the discovery of neutron rich light nuclei, 
it was found that the BW is severely inadequate for light mass nuclei away from the valley of stability. From a 
phenomenological search a more generalized BW formula was predicted by Samanta and Adhikari [4] in which the 
binding energy is given by, 

B(A, Z) = a v A - a s A 2 ' 3 - a c Z(Z - 1)/A^ 3 - a sym (A - 2Z) 2 /[(1 + e- A ' k )A] + S new , (4) 

where, 

8 new = (1 - e- A ' c % c = 30, k = 17, (5) 

and the constant a v = 15.777MeV [7]. The other constants remain the same and 6 is the old pairing term. 

We made five parameter searches for the BW mass formula using Audi-Wapstra 2003 [8] mass table by minimizing 
X 2 /-/V and the root mean square deviation a for atomic mass excesses. The results are listed in Table-1 below. We 
find that both and x 2 /N and the root mean square deviation a, which are defined later in the text, can not be 
minimized simultaneously with the same set of parameters. These five parameters obtained from x 2 /N minimizations 
are the same as those used in BWM [7]. Once the five parameters are fixed from the y 2 '/N minimization, the rest two 
parameters c and k were extracted by matching the experimental binding energy versus neutron number curves from 
Li to Bi. For all the 2518 nuclei available in Audi-Wapstra 2003 [8] mass table from Li to Bi, excluding 12 C and 20 < 21 F 
for which the experimental mass excesses are either zero or very small, the values obtained are (x 2 /N)bw — 3.250, 
(x 2 /N)bwm = 1-434, (<j)bw = 4.711MeV and (<j) B wm = 3.788MeV. 

As BWM mass fomula does not have shell effect incorporated, it can not predict the masses of the magic nuclei and 
nuclei with deformities. Also, as there is no Wigncr term, it fails to predict the extra stability of the N=Z nuclei. To 
predict the masses of all the nuclei more accurately, we introduce the isotonic shift [5] in the mass calculations with 
BWM. The isotonic shift corrected binding energy, BE cor (A, Z), is defined as, 

BE cor (A, Z) = BE Th {A, Z)-A n , (6) 

where 

A„ = {\/n)Y^BETh{Ai, Zi) - BE expt (A u Z { % (7) 

i 

and n=3 is the number of points taken before the right or left end of the curve. The BEth (A, Z) is the theoretical 
binding energy. To achieve the most accurate treatment of the masses, we employ the recently available masses given 
in the Audi-Wapstra mass table, 2003 [8]. There still exists a large number of unknown masses for some nuclei, most 
of which have extreme proton to neutron ratios. For example 61 As is the heaviest and 40 Mg is the lightest isotones for 
N=28 listed in the Audi-Wapstra table 2003 [8]. To predict the binding energy for N=28 isotone heavier than 61 As, 
viz. 62 Se, we use the masses of 59 Ga(Z = 31), m Ge{Z = 32) and 61 As(Z = 33) and cqns. (6) and (7). Similarly, we 
compute A„ from the masses of 40 Mg(Z = 12), 41 Al(Z = 13), 42 Si(Z = 14) to predict the mass of 39 Na which is a 
N=28 isotone lighter than 40 Mg. To accomplish this computationally we find that the formula 

N s = 0.968051Z + 0.00658803Z 2 , (8) 

provides a reasonable description of the neutron number N s of a nucleus with atomic number Z on the line of beta- 
stability. For a particular (Z,N) nucleus, N s decides which side of the line of beta stability the chosen nucleus lies. 



TABLE I. The parameters for the Bethe-Weizsacker mass formula found by fitting the available atomic mass excesses of 
Audi-Wapstra 2003 [8] mass table. 



a v 


a s 


a c 


&sym 




Minimized 


MeV 


MeV 


MeV 


MeV 


MeV 


quantity 


15.777(53) 


18.34(17) 


.710(3) 


23.21(10) 


12(2) 


x7-/V=3-34 


15.258(20) 


16.26(6) 


.689(1) 


22.20(5) 


10.08(85) 


ct 2 =15.12 
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The one proton and two proton separation energies S\ p and S2 P defined as the energies required to remove one 
proton and two protons from a nucleus are given by, 

Sip = B(A, Z) — B(A — 1,Z — 1) and S 2p = B(A, Z) - B(A - 2, Z - 2) (9) 

respectively, while the one neutron and two neutron separation energies Si„ and S2n defined as the energies required 
to remove one neutron and two neutrons from a nucleus arc given by, 

Si„ = B(A, Z) — B(A — 1, Z) and S 2n = B(A, Z) - B(A - 2, Z) (fO) 

respectively. 

The atomic mass excess AM^z for a nucleus with mass number A and atomic number Z is given by 

AM AiZ = ZAm H + (A — Z)Am n - a el Z 239 - B(A, Z) (11) 

where Amn = m p + m e — u = 7.28898483 MeV and Am n — m n — u = 8.07131710 MeV, m pi m n and m e are the 
masses of proton, neutron and electron respectively, the atomic mass unit u is 1/12 the mass of 12 C atom, AMa,z 
is the atomic mass excess of an atom of mass number A and atomic number Z and the electronic binding energy 
constant a d = 1.433 x 10~ 5 MeV. 

In the present work the one neutron and one proton separation energies have been calculated using the newly 
extended Bethe-Weizsacker (BWM) mass formula given by eqn. (4) with the isotonic shift, called BWMIS. The last 
stable nucleus from which the removal of a single neutron (and any more) makes the one proton separation energy 
negative defines a proton drip line nucleus and the last stable nucleus to which addition of a single neutron (and any 
more) makes the one neutron separation energy negative defines a neutron drip line nucleus. Proton and neutron drip 
line nuclei, thus obtained, over the entire (N,Z) region define the predicted neutron and proton drip lines. 

The two quantities of importance that can be used for a relative comparison between different theoretical models 
are X 2 /N and the root mean square deviation a which are defined as 

X 2 /N = (1/A0 ^[(AM Tk - AM Ex )/AM Ex } 2 (12) 

and 

a 2 = (l/N)Y^[AM Th - AM Ex ] 2 (13) 

respectively, where AMrh is the theoretically calculated atomic mass excess and the summation extends to TV data 
points for which experimental atomic mass excesses are known. For the calculations of x 2 /iV and the root mean 
square deviation cr, experimental masses [8] of 3175 nuclei have been used for the cases of improved liquid drop 
model (ILDM) [5] and BWM. For the case of Myers-Swiatecki (MS) [9], out of 3175 experimental masses, theoretical 
masses of 23 nuclei are not available. Therefore, we use only experimental masses of 3152 nuclei for calculating the 
X 2 /N and a for MS masses. The values \ 2 / N obtained for the MS, BWM and ILDM are 1.398, 8.203 and 20.525 
respectively, while a for those are 0.865 MeV, 4.247 MeV and 3.061 MeV respectively When isotonic shifts have 
been incorporated in the BWM and ILDM mass formulae, the \ 2 /N reduces to 7.758 and 10.950 respectively while 
a becomes 1.645 MeV and 1.529 MeV for BWM and ILDM respectively. In case of BWM the isotonic shifts have 
been incorporated beyond A=30 since the BWM [4] was designed to improve the predictions of low mass region. It 
is worthwhile to mention here that a is a measure of absolute error while the % 2 /N is a measure of relative error. 
In many cases, such as reaction Q value or separation energy and dripline calculations, absolute errors involved in 
atomic mass excesses are the quantities of concern and a plays more important role than y 2 /N . On the contrary, 
X 2 /N signifies the relative, that is, uniform percentage error involved in the mass predictions. In fact, in BWMIS the 
value of x 2 /N reduces from 7.758 to 0.907 if just 2 nuclei, namely 20 F and 21 F, for which experimental mass excesses 
are of the order of 10~ 2 MeV, are excluded from the calculations. This happens as in the computation of the X 2 /N, 
the experimental mass excesses AM Ex appear in the denominator and any small deviation of AMth from AM Ex is 
blown up excessively when divided by a number close to zero. But it is pertinent to note that the exclusion of these 
two nuclei causes insignificant change in a from 1.645 MeV to 1.644 MeV. Moreover, exclusion of these two nuclei 
does not alter any of the parameter values listed in Table-1. However, the final calculations have been performed 
including all the nuclei. 

In the Table-2 below, using the BWM mass formula with isotonic shifts beyond A=30, the atomic number Z and 
the neutron number N of the neutron and proton drip line nuclei have been tabulated along with the corresponding 
values of their one proton and one neutron separation energies (with errors in the parentheses). Uncertainty of 0.1 
MeV in the calculation of binding energies have been used in the low mass region upto A=30 where isotonic shifts 
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have not been incorporated. The new calculations predict 31 F and s He to be stable as can be seen in the Table-2. 
The values given in Table-2 could be useful for planning experiments on mass measurements. In order to demonstrate 
the change of sign of the one proton and one neutron separation energies, the atomic number Z and the neutron 
number N of nuclei with one neutron less than the proton drip line nuclei and nuclei with one neutron beyond the 
neutron drip line nuclei have also been listed along with the corresponding values of their one proton and one neutron 
separation energies. 

TABLE II. Proton and neutron separation energies on and just beyond the proton and the neutron drip lines using the 
newly modified BW mass formula with isotonic shift. 

Proton S p S„ One S p S n Neutron S p S n One S p S n 

beyond beyond 
proton neutron 
Drip drip Drip drip 



Z,N MeV MeV Z,N MeV MeV Z,N MeV MeV Z,N MeV MeV 
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107.137 .52( 13) 9.53( 11) 

108.138 1.85( 13) 11.13( 11) 
109,155 .11( 13) 7.63( 10) 
110,141 .97( 13) 9.57( 11) 
111,160 .21( 13) 8.37( 10) 
112,144 .89( 13) 11.04( 10) 
113,165 .41 ( 13) 6.70( 10) 
114,148 .95( 13) 10.86( 10) 
115,170 .18( 13) 8.41( 10) 
116,150 .10( 13) 10.97( 10) 
117,174 .03( 13) 8.45( 10) 
118,154 .23( 13) 10.80( 10) 



107.136 -3.34( 13) 11.23( 11) 

108.137 -2.00( 13) 9.75( 11) 
109,154 -.44( 13) 8.65( 11) 
110,140 -2.18( 13) 11.24( 11) 
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112,143 -2.36( 13) 9.69( 10) 
113,164 -.13( 13) 7.90( 10) 
114,147 -2.35( 13) 9.52( 10) 
115,169 -.01( 13) 6.65( 10) 
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117,173 -.14( 13) 6.95( 10) 
118,153 -.27( 13) 9.48( 10) 
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112,282 18.69( 14) .04( 10) 

113,284 17.30( 14) .07( 10) 
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In Fig. 1, using the BWM mass formula, the atomic number Z has been plotted against the neutron number N for 
the proton and neutron dripline nuclei. The lines joining the proton dripline nuclei and the neutron dripline nuclei 
have been shown as the proton dripline and the neutron dripline respectively in the figure. The nuclei stable against 
decay via 0^, a or heavy particle emissions or electron capture or spontaneous fission [10] have been shown with 
small solid squares. There is no nucleus beyond the atomic number Z=83 which is stable against the decays via 0^, 
a, heavy particle emissions and electron capture and spontaneous fission. The hollow circles depict the proton and 
the neutron dripline nuclei using the BWM mass formula with isotonic mass shift corrections (BWMIS) beyond A=30 
while the solid triangles represent dripline nuclei using MS masses. In Fig. 2, the driplines have been calculated using 
the MS masses with the isotonic shifts (MSIS). The solid triangles represent dripline nuclei using MS masses while 
the hollow circles represent dripline nuclei using BWM mass formula with the isotonic shifts beyond A=30. Fig. 2 
shows that both the driplines with MS and MSIS are almost the same except at few scattered points. Some of the MS 
masses near both the driplines are not available. Specifically, MS-prcdictions for some even Z proton dripline nuclei 
are not available around N=90 and Z=75. As there is no simple formula by which the missing masses of MS can be 
calculated easily, they are missing in the figures. Interestingly, the driplines predicted by BWMIS are almost similar 
to those obtained by MSIS except at large N or Z values where we find that BWMIS gives better mass predictions, 
that is, closer to available experimental masses than MS. 

In summary, the atomic mass excesses have been calculated using the modified Bethe-Weizsacker mass formula 
with isotonic mass shifts incorporated and compared with the predictions of the mass formula of Myers- Swiatecki. 
It is found that the agreements between the experimental masses and the calculated ones are greatly improved by 
introducing the isotonic mass shifts. The x 2 /N= 0.907 predicted by BWMIS is comparable to the x 2 /N= 1.398 
predicted by MS. Neutron and proton separation energies have been calculated using the modified Bethe-Weizsacker 
mass formula with isotonic mass shifts. The last bound isotopes of neutron and proton rich nuclei have been identified 
and the neutron and proton driplines have been provided. The significant changes by incorporating isotonic mass shifts 
have been observed for the proton dripline near heavier mass region. It is pertinent to note that BWM is a simple 
mass formula compared to the mass formula of Myers and Swiatecki. Nevertheless, with isotonic shift incorporated the 
modified Bethe-Weizsacker mass formula not only predicts almost the same driplines as the MS formula but provides 
better mass predictions for some nuclei with higher N, Z values as well. These observations are important as they 
suggest that the simple BWMIS prescription would be useful for theoretical computations of masses which are used 
in various computational codes. 
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FIG. 1. The atomic number Z has been plotted against the neutron number N as solid squares for the stable nuclei and the 
driplines have been calculated using the BWM mass formula. The solid triangles represent dripline nuclei using MS masses 
while the hollow circles represent dripline nuclei using BWM mass formula with the isotonic shifts beyond A=30. 
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FIG. 2. The driplines have been calculated using the MS masses with the isotonic shifts. The solid triangles represent 
dripline nuclei using MS masses while the hollow circles represent dripline nuclei using BWM mass formula with the isotonic 
shifts beyond A=30. 
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